Statistics of wave interactions in nonlinear disordered systems 
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We study the properties of mode-mode interactions for waves propagating in nonlinear disordered 
one-dimensional systems. We focus on i) the localization volume of a mode which defines the number 
of interacting partner modes, ii) the overlap integrals which determine the interaction strength, iii) 
the average spacing between eigenvalues of interacting modes, which sets a scale for the nonhnearity 
strength, and iv) resonance probabilities of interacting modes. Our results are discussed in the light 
of recent studies on spreading of wave packets in disordered nonlinear systems, and are related to 
the quantum many body problem in a random chain. 
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I. INTRODUCTION 

In the absence of nonlinearity (or many-body in- 
teractions in quantum systems) all eigenstates in one- 
dimensional random lattices with disorder are spatially 
localized. This is Anderson localization flj, which has 
been discovered fifty years ago in disordered crystals as 
a localization of electronic wavefunction. It can be in- 
terpreted as an interference effect between multiple scat- 
terings of the electron on random defects of the poten- 
tial. Recent experiments on the observation of Ander- 
son localization were performed with light propagation 
in spatially random optical media [2, i^] , with noninter- 
acting Bose-Einstein condensates expanding in random 
optical potentials jj, la] , and with wave localization in a 
microwave cavity filled with randomly distributed scat- 
terers |6J. 

In many situations nonlinear terms in the wave equa- 
tions (respectively, many body interaction terms in quan- 
tum systems) have to be included. Thus, a fundamental 
question which has attracted the attention of many re- 
searchers is what happens to an initial excitation of arbi- 
trary shape in a nonlinear disordered lattice. Nonlinear- 
ity renormalizes excitation frequencies, thereby inducing 
interaction between NMs. Numerical studies show that 
wave packets spread subdiffusively and Anderson local- 
ization is destroyed |7HlO| . In the regime of strong non- 
linearity, far from where it can be treated perturbatively, 
new localization effects of selftrapping occur [11| . A the- 
oretical explanation of the subdiffusive spreading was of- 
fered in Refs. [3, |^, [i2|. It is based on the fact that the 
considered models are in general nonintegrable. There- 
fore deterministic chaos will lead to an incoherent spread- 
ing. Estimates of the excitation transfer rate across the 
packet tail are obtained by calculating probabilities of 
mode-mode resonances inside the packet. Some predic- 
tions of this approach include the effect of different de- 
grees of nonlinearity and were successfully tested in [13| . 

In this work we study the statistical properties of 
mode-mode interactions. We focus on i) the localization 
volume of a mode which defines the number of interacting 




FIG. 1: (Color online) Localization length ^ versus normalized 
eigenvalue for W=0.5, 1, 2, 4 (from top to bottom). Inset: 
zoom for W=0.5 around the bandwidth center. 



partner modes, ii) the overlap integrals which determine 
the interaction strength, iii) the average spacing between 
eigenvalues of interacting modes which sets a scale for 
the nonlinearity strength, and iv) resonance probabili- 
ties of interacting modes. We discuss the results in the 
light of recent studies [3-1^, [lll - [l3j on spreading of wave 
packets in disordered nonlinear systems, and relate our 
findings to the quantum two interacting particle problem 
in a random chain. 



II. NONLINEAR SCHRODINGER CHAIN 

We consider the disordered discrete nonlinear 
Schrodinger equation (DNLS) with the Hamiltonian 



'Hd=Y.^i\M'' 



/? 



\M''-{^i+iri+ri+M- (1) 



Here ipi are complex variables, / are lattice site indices 
and /3 > is the nonlinearity strength. The random 
on-site energies e/ are chosen uniformly from the interval 
[— ■^, ^], with W denoting the disorder strength. The 

equations of motion are generated by tpi = dHo/diitp*): 



iipi = ei^i + Plipil^ipi - -ipi+i - Tpi-i 



(2) 



Eq. ([2]) conserves the energy ([T]) and the norm S = 
^i IV'iP- Varying the norm of an initial wave packet 
is strictly equivalent to varying /3, therefore we choose 
5 = 1. Note that Eq. ([2|) is used to qualitatively de- 
scribe the evolution of a dilute Bose-Einstein condensate 
trapped into a deep periodic potential [J], and also the 
evolution of a light wave in disordered one-dimensional 
waveguide lattices with cubic Kerr-type nonlinearity (un- 
der the paraxial approximation) [3j . 

For /? = 0, Eq. ([T]) with tpi = Ai exp(— iAi) reduces to 
the eigenvalue problem 



XAi = eiAi - {Ai+i + Ai^i 



(3) 



The width of the eigenfrequency spectrum Xi, of ([3]) is 
A = W+4: with X^ e [-2-^,2 + ^]. The normahzed 
eigenvectors A,yj {J2i^ti ~ 1) ^^"^ ^^e NMs, and the 
eigenvalues X^, are the frequencies of the NMs. We order 
the NMs in space by increasing value of the center-of- 
norm coordinate X^, = J^i^^t i- 

The equations of motion of ([T|) in normal mode space 
read 



itj)^ = X„cj)^ +;9 y^ /, 



U,l'i,U-2.l'3'rviVl'2^V3 



with the overlap integrals 



(4) 



(5) 



The variables <f>i, determine the complex time-dependent 
amplitudes of the NMs. 



III. PROPERTIES OF NORMAL MODES 



A. Localization length, volume and participation 
number 



The asymptotic spatial decay of an eigenvector is given 
by A^,i ~ e"'/^'' , where ^i^ is the localization length of a 
mode V with the eigenvalue Ai^. We calculate the average 
^^ at a given energy using the standard transfer matrix 
approach [14] and show the results in Fig.[TJ As expected, 
the most extended modes correspond to the bandwidth 
center with C(A = 0, W^) « lOO/PF^ for VF < 4 [H (see 
Fig. [2]) . In what follows we refer only to the localization 
length near the bandwidth center. We also observe a 
small peak at A = for W <1. The smaller the disorder 




FIG. 2: (Color online) Average localization volume V [(r)ed], 
participation number P [(g)reen)] and localization length ^ 
[(b)lue] of NMs with eigenvalues near the bandwidth center 
versus strength of disorder W . Dashed lines are estimated 
asymptotics for V and ^ at small disorder strengths, respec- 
tively, 330/iy^ and 100/1^^ 



for W = 0.5, the additional peak height is about 8% of 
the total value. Its magnitude will not exceed roughly 
10% for W — >• 0, as first discussed in Refs. [l^- The ori- 
gin of this anomaly is the deviation from single parame- 
ter scaling due to the symmetry ^;(A) ~ (— l)'A/(— A) at 
W = Q ISj. 

Next we estimate the number of NMs which interact 
with a chosen mode v. This is assumed to be equivalent 
to estimating the number of sites, where the norm density 
\A^^i p of the eigenvector is not exponentially small. This 
number is coined localization volume, Vy. Vi, is related to 
the localization length ^^^ , though quantitatively the two 
quantities might differ. We consider two ways of estimat- 
ing this number. A widely used quantity is the participa- 
tion number p^ — 1/ ^^ A'^ ^ . It is a measure of the inho- 
mogeneity of the distribution of eigenvector amplitudes 
in real space. Another quantity is the effective distance 
between the exponential tails of the eigenvector, which is 

|2 



given by yl2m 



,(>-) 



+ 1, where ?77,2 



Eiix.-ir\A,A 



is the second moment of the norm density distribution. 
Both quantities yield the exact width of a flat and com- 
pactly distributed norm density distribution. However 
when fluctuations are included, p^ will be reduced and 
underestimates the correct volume, while the effective 



12m 



(!') 



1 



strength, the more pronounced the peak is. For instance 



distance does not. Therefore, we use V^, - 
as a measure of the localization volume. 

We calculate numerically the average localization vol- 
ume V = Vi, and participation number p = p^, of NMs as 
a function of W (see Fig. [5]). For this purpose we fix the 
strength of disorder W and take a chain which is much 
longer than £,{W). We calculate V and p, taking into 
account only those modes whose eigenvalues are located 
near the bandwidth center and which are not close to 
the boundaries of a lattice. Then, we take another re- 



alization and repeat the procedure. Finally, we perform 
the averaging with respect to different disorder realiza- 
tions. We find that the localization volume V scales on 
average as 3.3i^ for weak disorder (see the dashed lines 
on Fig. [2]), and tends to y = 1 in the limit of strong 
disorder. We also note that the participation number is 
almost identical to the localization length for weak dis- 
order, and therefore misses the localization volume by a 
factor of three. 

Our numerical results indicate that multi-humped 
NMs are rare and do not affect the statistical results on 
V. These multi-humped NMs are localized over a num- 
ber of lattice sites which can be at far distance from each 
other in real space (see e.g. [131 )• ^^ they were statis- 
tically relevant, the second moments would be overes- 
timated. However, we find that in the limit of strong 
disorder W —?' oo the localization volume V —> 1, there- 
fore multi-humped states do not significantly contribute 
(although we admit that there might be a measurable 
contribution in this limit of strong disorder) . In the limit 
of weak disorder we observe that V, P, £, scale in the same 
way with the disorder strength, therefore we can exclude 
any significant statistical relevance of multi-humped NMs 
in this regime. 



B. Overlap integrals 

We study statistical properties of the absolute values 
of the overlap integrals ([5]) perturbatively for weak disor- 
der, and numerically using two different methods. In par- 
ticular we aim at estimating the average absolute value 
of these overlap integrals for NMs which are interact- 
ing with each other within the range of one localization 
volume, in order to exclude statistically irrelevant expo- 
nentially weak interactions of distant NMs. Note, that in 
following notations, the absolute value is omitted for the 
sake of simplicity. In order to avoid multiple repetitions 
in ([5]) we use I'a > t^2 > J^i > ^• 

Perturbative calculations. Let us consider a chain with 
finite size N and fixed boundary conditions: 



XAi = WiiAi - Ai+i - Ai^i 



(6) 



I = 1,...,N, -00 = ipN+i = and e are random un- 
correlated numbers evenly distributed over the interval 
[—1/2, +1/2]. For W — the canonical transformation 
to standing waves 



A,^ 



N 



N ■ 



tEq 



q-^ql ; ^ql 



sm 



irql 

N + 1 



(7) 



yields eigenvectors Ag^ = \ j^riSqi- Eq. ^ transforms 
to 



N 



ACj/q — ^q^q I /t^ / -ti-pq(oJp , /^ — 



2W 



p=i 



VTTTT 



, (8) 



with the coupling 



K„ 



Vn~ 



N 

= ^eiSgiSp 
1=1 



(9) 



which mixes standing waves with the eigenvalues Xq = 
2 cos {TTq/{N -1-1)) in the presence of disorder. 
The overlap integral 



/, 



qi,q2-q3-q4. 



N 

1=1 



lAq2,lAq.^lAq^l 



(10) 



at VF = will be zero for all combinations of indices 
except if a selection rule is satisfied ^M,- It is enough to 
replace this rule by 94 = ±qi ± (72 ± 93- In short we will 
denote by q^ a mode number which satisfies the selection 
rule for a given triplet of mode numbers {qi,q2, 93)- The 
selection rule applies to N^ overlap integrals /g ~ 1/A^. 
The other iV^ integrals /i = 0. Therefore, the average 
overlap integral becomes {I){W = 0) ~ 1/N^. 

Let us estimate the corrections to this average when 
disorder is added. We first consider integrals /i which 
were strictly zero at the limit W = 0. We perform a 
perturbation calculation {W small) for a mode (74 such 



that Qq = Qg 

Straightforward calculation gives (see also il9j 



with QfJ = 1 and g^X^ - 0. 






K. 



9,94 



A„ 



A„ 



(11) 



Assuming now a triplet of modes (gi, 52, 93) is given, and 
that 54 ^ 94, the first order nonzero correction to the 
corresponding overlap integral reads 



/, 



91,92, 93 -94 



(N 



N 



Sqxl^q^l^qal 7 , -""9 

1=1 q-^qi 



9ql 



'^q 

fl2) 



We started with q^ ^ 54, but in the presence of disor- 
der the mode with number q = q^ will become excited. 
Therefore, after summation over I in (|12p we find 



K, 



'91,92,93,94 



[N + 1) Xq, - Xq 



(13) 



Note that the indices (91,92,93) are implictely hidden in 
the quantity 94. In order to estimate the average, we 
have to take the absolute value of ()13|) . to sum over each 
index 9^, i = 1,2,3,4 and each time to divide by N. 
Let us perform the averaging over 94. The denominator 
Xq, — Xq, will become of the order of 1/iV when 94 is close 
to 94. Replacing the sum by an integral, we estimate 



1 



/ , l-'9l, 92, 93,94! 

94 #94 



iV+1 



In(iV)- 



1^^9-4.9-4! 



i(^9"4/iV)| 



(14) 



Since the disorder average (Kpq) — and its variance is 
finite (i.e. not depending on N) the final averaging over 

qi,q2,q3 yields 



(/i) r^ Kln{N)/N. 



(15) 



The overlap integrals Iq for 94 = 54 were of the order of 
1/A^ for W — 0. It is straightforward to obtain that the 
disorder induced correction will be of the same order as 
P^ . which is still smaller than the unperturbed value. 

Thus the average value of (/) up to the first order 
of perturbation in W is given by (remember that k = 

2W \ 

^/W+l' 



(/) - l/iV^ ^ aWln{N)/N^/^, 



(16) 



where a is some constant independent on the system's pa- 
rameters. One can conclude from Eq. ([T5)) . that for small 
enough W such that W < {\/N \n{N))~^ , the first term 
prevails and the total average integral is (/) ~ 1/iV^. In 
the opposite case, when W > {y/Nln{N))^^ , the sec- 
ond term in Eq. (J16p dominates and, as a result, we get 
(/) ~ W/7V-3/2ln(Ar). 



Note that the perturbed eigenvectors A 



A 



(0) 



kA'-^i' + ..., given by A'^/ = Y,^^^ ^9,p a^^ do not yield 
logarithmic divergence, since - at variance to the overlap 
integrals - no absolute values are taken, and the two log- 
arithms obtained from integrating to the left and right 
of q are cancelling each other due to opposite signs. 

For W —^ CO the NM eigenvectors become single site 
profile, and the overlap integrals tend to zero. There- 
fore, for a given size TV, the average overlap integral will 
start to increase with W for W > {^/N \n{N))~^ , reach a 
maximum at Wmax , and decay down to zero for infinitely 
strong disorder. It is reasonable to assume that the lo- 
calization volume ^{Wmax) ~ A^. In that case for small 
values of W we obtain 



{I) 



\n{V) 



-W^ln{W) 



(17) 



which is an estimate of the interaction strength of NMs 
within the spatial range of one localization volume. 

Numerical calculations. Method I. We fix a chain size 
N and calculate the average value of the overlap integrals 
taking into account all integrals ([5]). Then, we take an- 
other realization and repeat the procedure. Finally, we 
perform the averaging with respect to different disorder 
realizations. Each averaged integral (/) is a function of 
W. As derived in the above perturbation approach, it 
has a maximum value {I){Wmax) at a certain Wmax, as 
is shown in the inset of Fig. [3] for A^ = 40. Now we vary 
the chain size N, and repeat the procedure. In Fig. [3] we 
plot the maximum values of {I){Wmax) as a function of 
Wmax (red curve). We find that £,{Wmax)lN « 8/3, as 
expected in the above perturbation approach. We also 
find that for large N the data can be fitted with the 
power law (/) « 0.0034W, with a = 3.40 ± 0.02 (see 
Fig. [3]). The fit was done using different number of nu- 
merical points (from 3 to 14) starting from the small- 
est W. In all cases the RMS relative error was better 
than lO^'^. We expect that this method will overesti- 
mate the corresponding prefactors. This is due to the 
fact that N « Q.i^£^{Wmax) and therefore states overlap 




FIG. 3: (Color online) Average integrals (J) versus strength 
of disorder W using method I [(r)ed] and method II [(b)lue]. 
Dashed fitting lines are: 3.4- 10"^VF^* (upper line) and 3.84- 
10-5M/3.4 (lo.^^gj. ijjjg) ij^gg^. i^i^ versus W for A^ = 40. For 

the averaging, 400 disorder realizations were used. 



more strongly than in an extended system, as seen in the 
next method. 

Method 11. We fix the strength of disorder W and 
choose a chain size N ^ ^{W). We select a middle part 
of a smaller size (core) of the width L, and do not con- 
sider the edges in order to avoid boundary effects. We 
use only modes within the core {v = 1..L). For each 
mode V we calculate its localization volume Vi, . Now we 
consider only NMs which happen to reside in a corre- 
sponding neighborhood, i.e. we select V,y/2 modes from 
the right and left (in case V^, is odd, one mode is ran- 
domly taken from left or right in addition). Therefore, 
we have defined a subset of NMs which interact with the 
i/th NM. We calculate all overlap integrals for this sub- 
set. Then, we move on to the next reference mode from 
the core. This procedure is performed for all NMs from 
the core, for many realizations. For small W the data 
can be fitted with the power law (/) « 3.84 • lO-'^M^^'' 
(see Fig. [3|) . Note that both methods yield the same 
exponents. Note also that we lack more data to distin- 
guish between the numerically found law W^"^ and the 
perturbation result — VF*ln(l/F). 



IV. FREQUENCY SCALES 

There are two frequency scales set by the linear equa- 
tions ([3]) : the average spacing d of NMs within the range 
of a localization volume and the width of the spectrum 
^ 0jl9j 12]. The two scales d < A determine the packet 
evolution details in the presence of nonlinearity. In or- 
der to calculate the average spacing and its distributions 
numerically, we fix the strength of disorder W and take 
a chain which is much longer than y(iy). We select a 
middle part of a smaller size (core) , and do not consider 




FIG. 4: (Color online) PDFs of eigenvalue spacings s for 
W=4,10,20 [(r)ed,(g)reen,(b)lue]. (O)range curve: PDF for 
W=4 and a short chain with N=10. Dashed curve: Wigner- 
Dyson distribution with the average spacing d ~ 0.59 [see 
Eq. dUl)]. 
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FIG. 5: (Color online) Orange (light gray) curve: numerically 
calculated average spacing d versus strength of disorder W 
(see Sec. IIV[) . Red (dark gray) curve: width of the spectrum 
A. Dashed curve: the fit of d in the limit of H^ — ^ by 
d = W^ /37. Dashed-doted curve: the fit of d in the limit of 
W -^oohy d = W/3 - 2. 



Typical probability density functions (PDF) of the 
spacings s are shown in Fig. U) For strong disorder 
W ^ 1 the relative contribution of small spacings to 
the PDF becomes smaller. The reason is that the local- 
ization volume tends to one and for each reference mode 
z/ we take into account only a single neighboring mode. 
As a result, spacings between eigenvalues increase. 

We also note that the computed PDFs are far from 
following a Wigner-Dyson distribution 



■KS 



W^5i-e-4<P 



(18) 



Especially for small W, such a distribution could be ex- 
pected due to large localization lengths. However we 
find systematic deviations towards a Poisson distribution 
with an enhancement of the probability density at small 
spacings. This is due to the fact that NMs overlap in 
general only partially in real space. The Wigner-Dyson 
distribution is recovered only for very short chains, when 
N < ^{W) [see Fig. |4]. In this case, all eigenmodes oc- 
cupy the same volume, and level repulsion is recovered 
as expected. Interestingly, similar level repulsion occurs 
in short resonators in one-dimensional random lasers. In 
that case the PDF of spacings between frequencies of the 
neighboring lasing modes tends also to the Wigner-Dyson 
distribution j2ll |. 

In Fig. [5] we plot the result for the average spacing d. 
In the limit of small W we estimate the average spacing 
as d = A/V oc W^. For large disorder strength W ^ 1 
the localization volume tends to one. Thus, only two 
modes form a subspace of each reference mode v (one 
of which is the mode i' itself). Therefore, the spacing 
can be calculated by considering two numbers (emulat- 
ing two corresponding eigenvalues) which are randomly 
distributed within the width of the spectrum A. The av- 
erage distance between these numbers a; and y, assuming 
that X > y is 



d = — / dx {x- y)dy, N 



dx / dy. (19) 
Jo 



It follows d = A/3. In Fig. [5] the two theoretical esti- 
mates are shown to be close to the numerical data. 



the edges. For each mode v within a core we form its 
subspace which consists of those modes which live in its 
localization volume Vi, [see Sec. IIII AIIII^ for details]. 
We take the eigenvalues of these modes (including the 
eigenvalue of v-ih mode), sort them and compute abso- 
lute values of spacings between them. Then, we proceed 
to the next reference mode from the core. This procedure 
is performed for all NMs from the core and for many re- 
alizations, such that we end up with a large number of 
spacings (usually of the order of 10^) [20| . 



V. RESONANCES 

When a wave packet spreads, its size grows, and the 
norm density inside the packet drops. Previous stud- 
ies indicate that this process of spreading is due to res- 
onances in mode-mode interactions. These resonances 
lead to chaotic dynamics inside the packet, and to a sub- 
sequent incoherent spreading. Let us estimate the num- 
ber of resonant modes in the packet. Excluding secular 
interactions, the amplitude of a NM with |0,y|^ — n^ is 
modified by a set of three other modes fl = (/ii,/i2,M3) 




FIG. 6: (Color online) Probability densities W(i?^,pf,) of res- 
onant NMs in linear-log scale (main figure) and log-log scale 
(inset). Disorder strength W = 4,6, 8, 10, 15 [(b)lack, (g)reen, 
(o)range, (v)iolet, (r)edl. 




FIG. 7: (Color online) The constant C as a function of W in 
linear-log scale when taking into account contributions from 
i) all combinations [(r)ed]; ii) quadruplets [(b)lue], and iii) 
triplets [(v)iolet]. 



in first order in /3 as ([4]) (see Q) 



PV' 



'^^l^/i2'^M3 






^iM 



\^ + '^Mi ~ '^P 



I^a 



(20) 
(21) 



The perturbation approach breaks down and resonances 
set in when y^nU < |(/)jy |. Since all considered NMs be- 
long to the packet, we assume their norms to be equal to 
n. 

We perform a statistical numerical analysis by com- 
puting the PDF of R^^p. For a given NM ly we ob- 
tain i?,y,/jo = min^i?;^^^. Collecting R,y,pg for many v 
and many disorder realizations, we find the probability 
density distribution 'W{Ri,_pg). We also analyze separate 
contributions from three different types of interactions, 
namely from quadruplets (all four modes are different). 



triplets (only three of four modes are different) and pairs 
(only two different modes participate in the interaction). 
For quadruplets all indices in (PT|) should be different, i.e. 
Ml i^v, ^2^^ V, ^3 ^v, ^i^ [12, Ml 7^ i^, Ml 7^ M3, M2 7^ 
/j,3. For triplets either ii\ — v such that 



(22) 





^,i^ ^^ 


A^2 2Aiy + A^3 


or [i2 - 


= /i3 with 






Rv,jl ^ 


\v — 2Ap2 + A^j 

J-iy,fJ.l,fJ.2,fJ.2 



(23) 



The remaining cases form the subset of pairs. 

The probability densities W{R,j^pa) of NMs being res- 
onant when taking into account all contributions are 
shown in Fig. [51 The main result is that W{R,y_pg — >■ 
0) -> C{W) ^ 0. The constant C drops with increas- 
ing disorder strength W (see Fig. [T]). We also calculate 
C by taking into account only quadruplets and triplets 
(see Fig. [7]). We find that for weak disorder the quadru- 
plet contributions are the dominant ones, while for strong 
disorder their contribution diminishes as compared to the 
triplet contribution. 

For small R the probability densities WiR) can be ap- 
proximated as 



WiR) « C(Ty)e-^(^)«. 



(24) 



The probability V for a mode, which is excited to a norm 
n, to be resonant at a given value of the interaction pa- 
rameter (3 is given by 



V 



pn 



W{R)dR « 1 - e 



'C0n 



(25) 



VI. DISCUSSION 

We have studied statistical properties of eigenvalues 
and eigenvectors of waves in disordered one-dimensional 
systems as a function of the disorder strength. We esti- 
mated the localization volume of a mode which defines 
the number of interacting partner modes. We obtained 
the dependence on the disorder strength of the overlap 
integrals which determine the interaction strength. We 
analyzed the statistics of level spacings of normal modes 
within one localization volume. Finally, we obtained dis- 
tribution functions for resonance probabilities of normal 
modes interacting in the presence of nonlinearity. Let us 
discuss some of the consequences of our findings. 



A. Overlap integrals 

In order to estimate the absolute value of the overlap 
integral ([5]) for modes within one localization volume for 
weak disorder, Shepelyansky 22] and Imry 23[ assumed 



that the sum extends roughly over the locahzation vol- equation in accordance with Q 



ume V, with each term in the sum A^iA^-^iA^^^iA^^i 
having a random sign. The absolute value of the eigen- 
vector is of the order of 1/V^^'^ due to normalization. 
Then ([5]) can be evaluated using the central limit theo- 
rem, for which the average absolute total value {I)rp ^ 
V~^l'^ r^ W^ . Our numerical finding (/) ~ V^3.4 dg^rly 
rules out the random sign resut W^ . As shown in the 
perturbation calculation in section III, the reason for the 
random sign failure is that NMs are similar to plane 
waves with definite phases on each lattice site (inside 
the localization volume). These phases enforce selection 
rules, which become strict in the very limit W = ^. 

While we can now exclude the random sign result 
W^ ^ we can not tell whether the numerical estimate 
(/) ~ Vl^3.4 jg correct, or the perturbation result (/) ^ 
— VK^ln(W^) will set in for small enough W . Ponomarev 
and Silvestrov 24] have also stressed the importance of 
phase correlations in Eq. ([5]) . A numerical calculation of 
the average of the squared overlap integral was performed 
by Frahm et al [H] for 1.4 < VF < 4 yielding (/) - VF^-^, 
in a good agreement with our numerical data. 

The random sign estimate W^ was taken to predict 
a strong increase of the localization length of two in- 
teracting particles in a one-dimensional random quan- 
tum chain [22, [23| . The two particle localization volume 
Vi , within a renormalization group approach, is given by 
V^jV ~ {Tf'V'^ where V is the single particle localization 
volume. For the random phase result, this yields V2 ~ V'^ 
[22 . I23I . l26j . We can clearly rule out such an outcome. 
Instead, we expect either V2 ^ V'^'^ (numerical data) or 
V2 ~ y In Y (perturbation approach) , which give a much 
weaker effect. These controversies call for more detailed 
investigations. 



B. Asymptotic spreading of wave packets in 
nonlinear chains 



According to a recent analysis of the spreading scenaria 
of wave packets [iJI , the only scale which separates differ- 
ent dynamical spreading regimes is the average spacing 
d. Therefore, the constant C from the previous section 
is inversely proportional to the mean level spacing: 



C ' 



(26) 



Following the theory developed in [3, [3, llJl for the 
asymptotic spereading, an exterior mode 0^ which is 
heated up by the packet obey the following evolution 



^0M « A^(/.^ + ^{I)V^V{^n)v?l'']{t), (27) 

where {f{t)f{t')) — S{t — t') ensures that f{t) has a con- 
tinuous frequency spectrum. Note, that here we also 
introduce the contribution of the overlap integrals es- 
timated as {I)V^ . Repeating the previous derivations 
0, M, [l3l ! we finally get the following expression for the 
asymptotic growth of the second moment of spreading 
wave packets in nolinear chains 



7712 



^4/3y8/3^j^2/3^1/3 



(28) 



From our numerical data for weak disorder it follows 
77^2 ~ ]4/-307^4/3^i/3^ while the perturbation approach 
yields 7712 - W-^'^{-\nWf/^l3'^/H^/^. The prefactor 
dependence of W is another intriguing test which awaits 
numerical verification. 



VII. CONCLUSION 

In conclusion, we performed a statistical analysis and 
calculated the average localization volume occupied by 
an eigenmode as a function of disorder strength which 
determines the average number of a nonexponentially in- 
teracting eigenmodes. Then, we calculated the frequency 
spacings of the normal modes which happen to interact in 
a nonexponentially weak way and their distributions and 
the average numerically. This result is very important 
for the classification of different regimes of wave packet 
spreading in the presence of nonlinearity. We also stud- 
ied statistical properties of the overlap integrals which 
determine the coupling strength between the interact- 
ing modes and, thus, influence properties of spreading. 
Finally, we estimated the number of resonant modes in 
the packet and proved that the most significant contri- 
bution to the spreading comes from the quadruplet and 
triplet resonances for small to moderate values of dis- 
order strengths, and from triplets for the case of large 
disorder. 
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